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We investigate entanglement dynamics in multipartite systems, establishing a quantitative concept 
of entanglement flow: both flow through individual particles, and flow along general networks of 
interacting particles. In the former case, the rate at which a particle can transmit entanglement 
is shown to depend on that particle's entanglement with the rest of the system. In the latter, we 
derive a set of entanglement rate equations, relating the rate of entanglement generation between 
two subsets of particles to the entanglement already present further back along the network. We 
use the rate equations to derive a lower bound on entanglement generation in qubit chains, and 
compare this to existing entanglement creation protocols. 
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I. INTRODUCTION 



New fields of physics often give rise to new physical 
quantities to study, and quantum information theory has 
proved a rich source of study material. As an amalgam of 
quantum mechanics and information theory, many of the 
new quantities are quantum analogues of familiar friends 
from classical information theory: the qubit, for instance, 
measures quantum information just as the bit measures 
classical information [lj. Other quantities have no ob- 
vious classical counterpart. The best-known example is 
entanglement. Originally seen as the phenomenon that 
epitomized quantum weirdness, it has become established 
over the last decade physical quantity, on a par 

with, say, energy. 

The analogy with energy can be pushed quite far: en- 
tanglement has a number of similar properties. Like en- 
ergy, entanglement can be quantified in a meaningful 
way 0], allowing us to say that one state is more en- 
tangled than another. Like energy, entanglement can be 
converted from one form to another [3( . And like energy, 
it is a resource that can be used to carry out useful tasks, 
such as teleportation Q. 

Until recently, work concentrated on understanding 
these static properties of entangled quantum states. Al- 
though we are some way from a complete understand- 
ing of entanglement statics, there has been significant 
progress: for instance, bipartite pure-state entanglement 
is now well understood. This begs the question: what 
happens if we allow the state to evolve? 

The move from entanglement statics to entanglement 
dynamics raises many new and interesting questions. 
How does entanglement evolve as particles interact 0? 
How good is aparticular interaction at creating entan- 
glement pJlllQ? More generally, how good is an inter- 
action at simulating various non-local processes 0, 0? 
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Or, turning this on its head, how 'non-local' is a given 
process (e.g. a quantum gate) 0,E1? This article ex- 
tends the first of these — how entanglement evolves as 
particles interact — to multipartite systems. 

The Schrodinger equation already implicitly describes 
the complete dynamics of a quantum system, but to gain 
insight into entanglement dynamics, we need equations 
that explicitly involve the entanglement of the system, 
without reference specific features of the Hamiltonian. 
One of the first steps along this path was taken by Diir 
et. al. who investigated the rate of entanglement gener- 
ation in two-qubit systems [f| . They derived an equation 
relating the rate of entanglement creation to the existing 
entanglement in the system, along with a factor depend- 
ing on the form and strength of the interaction. This 
latter led to a pleasingly simple quantity measuring the 
entanglement generating capacity of two-qubit interac- 
tions |E0- 

In a system of two particles coupled by a Hamiltonian, 
the only entanglement dynamics that can take place is 
creation of entanglement between the two particles. A 
simple tripartite system already raises other interesting 
questions. For instance, in a chain of three particles, 
how does entanglement 'flow' through the middle one? 
Surprisingly, we showed in previous work |12| that, in 
just such a chain, entanglement can be created between 
the two end particles, without the middle particle ever 
becoming entangled. This would seem to put an end to 
notions of entanglement 'flow'. However, we also gave a 
simple proof that this phenomenon is only possible for 
mixed initial states; for pure states, the middle particle 
necessarily becomes entangled during the evolution. 

This suggests there is a connection between pure-state 
entanglement of a mediating particle and entanglement 
flow through that particle: if it is not entangled, no en- 
tanglement flows. In section [HJ we show that there is 
indeed a quantitative relation describing how the entan- 
glement of a particle with the rest of the system limits 
the flow of entanglement through that particle. We first 
consider a three-qubit system, before dealing with gen- 
eral systems. The concept of entanglement flow through 
particles is therefore put on a quantitative footing for 
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systems in pure states. This contrasts strongly with the 
mixed-state case, in which entanglement can seemingly 
'tunnel' through mediating systems. 

Flow through individual particles is one aspect of en- 
tanglement dynamics in multipartite systems. But in a 
network of many interacting particles, we may also be in- 
terested in how entanglement flows along the whole net- 
work. We develop these ideas in the second half of this 
article. 

The inspiration is loosely based on the Arrhenius equa- 
tions for chemical reactions. The reaction mechanism of a 
chemical reaction describes the steps by which reactants 
are transformed, via successive intermediate compounds, 
into the final products. The rate at which a compound is 
produced depends on the amounts of its immediate pre- 
cursors that are present. Thus the complete reaction is 
described by a set of coupled rate equations, one for each 
step in the reaction mechanism. 

In section IIIII we derive a set of differential equations 
describing entanglement flow, analogous to the rate equa- 
tions for a chemical reaction. The rate at which entangle- 
ment is generated between two sets of particles is shown 
to depend on the amount of entanglement already present 
further back along the network. The entanglement dy- 
namics of the complete system is described by a coupled 
set of such entanglement rate equations, one for each step 
in the interaction network. 

Unlike the equations describing flow through individ- 
ual particles, these entanglement rate equations apply 
equally well to both pure and mixed states. Therefore, 
they establish a concept of entanglement flow along gen- 
eral networks of interacting particles (even though the 
concept of flow through individual particles in the sys- 
tem may be meaningless). 

In section ITVl we apply our new understanding of en- 
tanglement flow to investigate entanglement generation 
in chains of interacting particles. First, we briefly re- 
view some existing entanglement generation protocols for 
qubit chains, in the context of the rate equations derived 
in section UTTl Finally, we use the rate equations to prove 
a universal lower bound on the time it takes to create en- 
tanglement, or more precisely the scaling of this with the 
length of the chain (the results can easily be extended to 
general networks). 



II. FLOW THROUGH PARTICLES 

In this section, we will investigate entanglement flow 
through mediating particles. Specifically, we will con- 
sider flow through the middle particle in tripartite chains. 
The results of [lj show that this concept does not make 
sense if the whole system is in a mixed state. But for 
pure states, the rate at which entanglement is generated 
between the end particles is indeed zero if the middle 
particle is not entangled. 

The latter is suggestive: is there a general quantita- 
tive relationship between the entanglement of a particle, 



and the rate at which entanglement can flow through it, 
for systems in pure states? If the middle particle is only 
slightly entangled, does entanglement flow only slowly? 
We will derive just such a relationship, first for the sim- 
plest tripartite system: a thrcc-qubit chain, then for gen- 
eral tripartite chains. 

When investigating (bipartite) entanglement flow 
through mediating particles in more general settings, the 
system can always be described as a tripartite chain: the 
mediating particles form one party, and the sets of par- 
ticles each side, which are becoming entangled, form the 
other two. Thus the equation for tripartite chains can in 
fact be applied generally to describe entanglement flow 
through mediating systems. 

A. The three-qubit chain 

Consider a chain of three qubits, labeled a, b, and c, 
with nearest neighbour interactions described by Hamil- 
tonians H a i, and Hi, c . We will restrict the overall state of 
the system, \ip) abe to be pure. However, the reduced state 
of the two end qubits, p ac , need not remain pure during 
the evolution (if b is to become entangled at any point, 
p ac will necessarily become mixed). To quantify the en- 
tanglement between a and c, we need an entanglement 
measure valid for mixed states. The natural choice is 
the concurrence [l3l | . Though it is an entanglement mea- 
sure in its own right, its interest lies in its is equivalence 
to one of the important, physically meaningful entangle- 
ment measures: the entanglement of formation |l4j . 

We can write the overall state of the system in its 
Schmidt decomposition with respect to the partition 
(6 : ac): \ip) abc ^ = Ai \tpi) ac \ X i) b + A 2 |y? 2 ) oc |xa> 6 < The 
Schmidt coefficients, Ai and A2, determine the non-local 
properties of the state with respect to this partition, in- 
cluding entanglement of the middle qubit b with the rest. 
Meanwhile, the entanglement of the reduced state of the 
end two qubits, p ac , can be measured by the concur- 
rence, denoted C ac . Following 0, the state of par- 
ticles a and b can be represented by a 4 x 2 matrix 
X = (Ai \(pi) , X2 \<f2})- The concurrence can be calcu- 
lated from the singular values ?i > q of A — X T HX , 
where S = a y ® a y : C ac = ft — ft = y/trA^A — 2| det A\. 
Taking the time-derivative of this and simplifying the 
resulting exact expression (see Appendix El f° r details) 
leads to the following bound on the entanglement rate: 

AC 2 

<8||Jf||AiA 2 . 

at 

The factor \\H\\ — ||-ff a &||i + ||-£f&c||i measures the 
strengths of the interactions. = Y^ij l-^yl denotes 

the l\ norm, where the Hamiltonians are written in the 
product basis H = HijUi ® aj, and the Pauli ma- 
trices 0X2,3 — <^x,y,z are defined in the Schmidt basis 
{|Xi) |X2)}- (Local terms Ui ® 1 or 1 ® <jj in the Hamil- 
tonian can not alter the entanglement of the system, so 
do not contribute). 
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In the context of entanglement dynamics, the impor- 
tant part of the relation is the product of Schmidt coeffi- 
cients A1A2, which is a pure-state bipartite entanglement 
measure. (In fact, up to a numerical factor, it is the 
concurrence.) Thus the differential equation tells us that 
the entanglement of the middle qubit limits entanglement 
generation between the end qubits: not only must b be 
entangled for entanglement to be generated between a 
and c (precisely what was shown not to hold for mixed 
states in 0), but the rate at which it is generated can 
be larger the more entangled b is. 

At first sight, the inequality may appear too weak, as 
it does not seem to imply that the derivative is zero once 
qubits a and c are maximally entangled. However, C ac 
and A1A2 are not independent quantities. When a and c 
are maximally entangled, they can not be entangled with 
anything else, thus A1A2 = 0, and the derivative is zero 
after all. 

A complete quantitative description of entanglement 
creation in the three-qubit chain would require an equa- 
tion describing the evolution of the Schmidt coefficients. 
However, (b : ac) forms a bipartite, pure-state system. 
Entanglement creation in bipartite systems and the evo- 
lution of the Schmidt coefficients has been investigated 
in 0. 



B. Fidelities and entangled fractions 

The three-qubit result can not directly be extended 
to higher dimensional systems. Whilst we can restrict 
the system to pure states for the same reasons as in the 
three-qubit case, the reduced density matrix of the two 
end particles can again become mixed during the evo- 
lution. And no closed-form expression is known for the 
entanglement of formation of mixed states, other than in 
the two-qubit case. 

Before turning to higher-dimensional systems, it is in- 
structive to consider more carefully the setting in which 
we wish to investigate entanglement dynamics. Entangle- 
ment measures are defined in the LOCC paradigm: lo- 
cal operations and classical communication (LOCC) can 
only decrease the entanglement of a state. This is the nat- 
ural paradigm when thinking about entanglement from 
an information-theorist's point of view, in which entan- 
gled states are shared between different parties who are 
free to act locally on their part of the state. 

But we are considering entanglement dynamics from 
a physical standpoint. In a system of interacting parti- 
cles, it is not clear what classical communication means. 
Any transfer of classical information between particles 
would still have to take place via the (quantum) inter- 
actions. It could be argued that it makes more sense in 
this context to define entanglement in the local-unitary 
paradigm: any change to a state due to local terms in the 
interaction Hamiltonian should not change the entangle- 
ment. 

A physical way of measuring entanglement in this 



paradigm is to use the fidelity ^(|, which measures the 
distance between states |34j . The entangled fraction of 
a state p is then defined as the maximum fidelity with a 
maximally entangled state: 

F{p) ■= max (<P\p\<P), 
\<j>)eM.E. 

where the maximization is over all maximally entangled 
states \cf>} in the bipartite Hilbert space of p. (For two- 
qubit states, it is also called the singlet fraction.) It 
measures how close a given state is to any maximally 
entangled state, and is invariant under local unitary op- 
erations, as required. 

The entangled fraction is also an experimentally rele- 
vant quantity. When trying to engineer an evolution to 
produce a particular state (a highly entangled one, for 
instance), we want to know how close the actual state is 
to the desired one — precisely what is measured by the 
fidelity. For example, in teleportation experiments, it is 
the entangled fraction of the entangled pair that deter- 
mines how close the teleported state is to the original • 

Therefore, in the remainder of this article, we will 
consider evolution of the entangled fraction and related 
quantities. Though it is a well- motivated quantity to 
study in its own right, it can also be used to give upper 
and lower bounds on entanglement measures such as the 
concurrence 01 (and hence entanglement of formation) . 
In particular, if a state is separable, its entangled fraction 
is less than or equal to 1/n (with n the dimension of the 
smaller of the two Hilbert spaces making up the bipartite 
space). Whereas if (and only if) the entangled fraction 
is equal to one, the state must be maximally entangled. 

In the final section, we will use our results to derive 
bounds on how long it takes to entangle particles when 
the system starts in a separable state. In this context, 
any quantity that takes different values for separable and 
maximally entangled states is equally good in principle: 
we can bound the time required to change from one value 
to the other. The entangled fraction, for example, must 
increase from 1/n to 1 . 

C. General tripartite chains 

The tripartite chain is a prototype for all indirect (bi- 
partite) entanglement creation. We can always divide a 
system into three: two systems that are being entangled, 
and everything else lumped into one mediating system. 
We can then investigate entanglement flow through this 
mediating system. 

In a general tripartite chain, consisting of systems A, 
B and C of arbitrary dimension, interacting by nearest- 
neighbour interactions Hab and Hbc, the Schmidt de- 
composition has the form \4>)abc = Si ^ I^Mac I*)b> 
where we sort the Schmidt coefficients A^ in descending 
order. By re-expressing the entangled fraction as a max- 
imization over purifications using Uhlmann's Theorem 
(see Appendix 0), we can derive an exact expression for 
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the time derivative of the entangled fraction. Simplifying 
the exact result to separate out the entanglement depen- 
dence yields a relation analogous to the three-qubit case 
(see Appendix IU1 for details): 

F(pac) < 2 \H\ VH^c) (2 W 3 - A?) 

ij 

Again, the factor \H\ = \Hab\qo + \F[ab\ 00 measures the 
interaction strengths, independent of the system state 
(l-Hloo = maxjj \Hij \ denotes the l^ norm). 

The quantity in brackets is closely related to the en- 
tangled fraction of \tp) ABC in the (B : AC) partition: 
F{\tjj)(%p\) = ~ J2ij ^i^j (with n the smaller of the di- 
mensions of B and AC). Subtracting Xf re-scales this 
entangled fraction so that it is zero when the state is sep- 
arable. Therefore, the entanglement of B with the rest 
of the system limits the rate at which entanglement can 
flow through B. As in the three-qubit case, the deriva- 
tive implicitly goes to zero when systems A and C be- 
come maximally entangled, since they can not then be 
entangled with B. 



III. FLOW ALONG NETWORKS 

In the previous section, we examined entanglement 
flow through the middle particle in a tripartite chain, and 
noted that the results can be applied to flow through in- 
dividual particles in general systems, by viewing the sys- 
tem as a tripartite chain. However, in a large multipar- 
tite system, this approach means lumping many particles 
together into single composite particles, hiding much of 
the entanglement dynamics. Can we more fully describe 
entanglement flow in networks of interacting particles? 

In this section, we derive a set of differential equations 
describing the entanglement dynamics, analogous to the 
rate equations for a chemical reaction. These show that 
the rate at which entanglement is created between two 
sets of particles depends on the existing entanglement 
further back along the network. Intuitively, this can be 
interpreted as entanglement flowing through the network. 



A. Generalized singlet fraction 

As in the previous section, we must first address the 
problem of how to measure entanglement in large sys- 
tems. Even before that, we must decide what entan- 
glement to measure, since multipartite systems provide 
a plethora of possibilities. What questions are we in- 
terested in investigating using our putative equations? 
Perhaps the most natural goal, given a system of many 
interacting particles, is to entangle a particular pair of 
them: the end qubits in a chain, for example. We will 
take this as our motivation for again considering entan- 
gled fractions of the two particles. We will also need to 



define a new fidelity-based quantity to measure bipartite 
entanglement embedded in larger systems. 

First note that, since any maximally entangled state 
can be reached by acting with local unitaries on a partic- 
ular maximally entangled state, we can of course maxi- 
mize over unitaries rather than states in the definition of 
the entangled fraction: 

F(p) = max (0| U\ g> U 1 b p ab U a ®U b \<j>) 

We can equally well think of the unitaries as acting on p 
rather than on the entangled state \4>) ab - This suggests an 
alternative interpretation of the singlet fraction: as the 
maximum fidelity with a particular maximally entangled 
state (e.g. the singlet) that can be achieved by acting 
with local unitaries. 

Based on this interpretation, we define the generalized 
singlet fraction, a measure of two-qubit entanglement for 
bipartite systems of arbitrary dimension (it can be ex- 
tended in the obvious way to measure general bipartite 
entanglement jl^l: 

F(pab) = max (<j>\ tv /ab (U A ® U B Pab U\ ® U B ) \<j>) , 

Ua,Ub 

(1) 

where a and b are qubit systems embedded in A and B re- 
spectively, \4>) ab is the singlet state, and the notation tr/ ob 
indicates the partial trace over all systems other than a 
and b. It measures the maximum fidelity with the singlet 
achievable by local unitaries. 

Note that, in two-qubit systems, this generalized sin- 
glet fraction reduces to the usual singlet fraction. For 
any system, it takes values between and 1, and for sep- 
arable states it is less than or equal to 1/2. Also, from 
the definition, if A and B are subsystems of A' and B', 
so that pab = tx/ AB {pA'B'), then F(pab) < F{p A 'B')- 

B. Entanglement rate equations 

We are now ready to state our main result: a set of cou- 
pled differential equations describing entanglement flow 
in networks of interacting particles. For simplicity, we as- 
sume that among the set of interacting particles S, there 
are at least two qubits a and 6, the premise being that 
we arc interested in entangling these. (The results can 
easily be generalized: see [jjj and Conclusions.) Let A 
and B be disjoint subsets of S. The equations describe 
the rate at which the generalized singlet fraction of pab 
can increase. 

Define A' and B' to be the sets of particles directly 
connected by an interaction to at least one particle in A 
or B respectively (i.e. A' is the set of particles at most 
'one-away' from A, thus A C A'; see FigQJ. If A' and B 1 
are disjoint (as in Fig.QJ, then the time derivative of the 
generalized singlet fraction is bounded by 

F(pab) < 2 ||JT|| VF(pab)VF(pa>b>) - F(p AB ), (2a) 
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FIG. 1: A network of interacting particles, showing inter- 
actions and sets denned in the entanglement rate equations. 
Interactions 'crossing the boundaries' of A or B are indicated 
by thicker lines. 

whilst if A' and B' have one or more particles in common, 
then 

F(pab) < 2 ||JT|| VF(pab)V1-F(pab)- (2b) 

The factor \\H\\ is a sum of strengths of those interactions 
that connect a particle in A or B to one outside A or B 
respectively (i.e. interactions that 'cross the boundary' 
of A or B; see Fig. 0: 

pf|| = £||i%IU, 

where || •||#<j denotes the Hilbert-Schmidt norm. 

The first step in the proof of these entanglement rate 
equations to rewrite the generalized singlet fraction Q in 
terms of purifications of pab using Uhlmann's theorem 
( Appendix iBjl . This leads to the following exact expres- 
sion for the derivative of the generalized singlet fraction: 

F(pab) = VF(pab)- max - £ (<p\ Hy \<p) 



where is a purification of pab, \x) is an extension 
of the singlet state to the Hilbert space of and 
\ip) =U\®U^ B |x). Using Lemma (Appendix |BJ, the 
terms inside the sum can be bounded by \ip) — 

h.c.) < 2 \\Hij || HS V(tr l^l) 2 - (tr(ReXy)) 2 , where 
= tr/jj |?/>)(c/?|. Finally, the quantities under 
the square-root can be related to generalized singlet 
fractions: (tr(ReXy)) 2 = F(pab) and tr(|Xy|) 2 < 
F(pA'B'), which concludes the proof. (The proof is given 
in full detail in Appendix [HJ . 

We can gain some insight into entanglement dynamics 
by considering the qualitative meaning of the rate equa- 
tions, before thinking about solving them. They divide 
a network of interacting particles into pairs of concentric 
sets, surrounding qubits a and b. For example, in Fig. ^ 
there are three such pairs: the qubits a and b themselves, 
the sets labeled A and B, and those labeled A' and B' . 
The rate equations tell us that entanglement must first 
build up between the largest sets, before it can cascade 
down successively smaller ones, finally reaching the two 



qubits (just as in a chemical reaction, intermediate com- 
pounds in the reaction mechanism must be created before 
the final product is reached). What is more, the rate at 
which the entanglement flows from one level to the next 
depends on the difference in entanglement between the 
two levels (somewhat like the rate of a reversible chem- 
ical reaction, which depends on the difference between 
the concentrations of reactants and products; or like flow 
in fluids, in which the flow rate depends on the pressure 
difference) . 

The number of pairs of sets is equal to half the 'inter- 
action distance' of the two qubits (rounded down to the 
nearest integer), i.e. half the smallest number of links in 
the network needed to connect a to b (in Fig. ^ their in- 
teraction distance is 5). A generalized singlet fraction can 
be defined on each pair of sets, along with an accompany- 
ing rate equation describing its evolution. Therefore any 
network has the same rate equations as a chain whose 
length is equal to the interaction distance, and whose 
interaction strengths along each link of the chain equal 
the factors \\H\\; all entanglement flow is equivalent to 
flow along a chain. This is qualitatively similar to results 
from quantum random walks, in which a quantum walk 
over a network is equivalent to a quantum walk along a 
chain 20]. 

Note that the factor \\H\ \ in the rate equations indis- 
criminately includes all interactions that cross the bound- 
ary. We might expect different interactions to contribute 
differently, depending on their location in the network. 
In fact, in Appendix [51 we derive a more general version 
of the rate equations, which accounts for each possible in- 
teraction pathway separately, and can therefore take into 
account the different roles different particles play in the 
entanglement dynamics, due to their differing connectiv- 
ity. The inequality in the corresponding rate equations 
is therefore tighter, but it leads to exponentially (in the 
number of particles) more equations describing the en- 
tanglement dynamics of a system, and gives a less intu- 
itive picture of entanglement flow. We will find that the 
simpler form given here is sufficient to derive a number 
of interesting results. 



C. Limits from the rate equations 

It is straightforward to prove inductively that the 
curves produced by saturating the inequalities in the rate 
equations I|2al2b|l constitute upper bounds on the evolu- 
tion of the generalized singlet fractions; i.e. the fastest 
possible evolution allowed by the rate equations is that 
which saturates the rate equations at each point in time. 

If the interaction distance between the qubits we intend 
to entangle is d, then the full set of rate equations involve 
\d/2\ generalized singlet fractions, which we will denote 
Fk(t), k = l...|_d/2_|. (|_»J denotes rounding down to 
the nearest integer.) We number them such that Fu/ 2 i 
is the singlet fraction of the two qubits. If we define 
Fa = 1, then the evolution of each Fk(t) is described by 
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equation JSaJ. 

Let fk(t) be the curves that saturate the rate equa- 
tions, i.e. fk{t) is the solution to 

fk = 2 \\H\\ \/Jk\/ fk-i — fk, 

with fo — 1. (For simplicity, we can take all coupling 
strengths \\H\\ to be 1.) Assume that fk(t) is an upper 
bound on F k (t), i.e. f k {t) > F k {t) for all t. If f k+1 (t) 
is not an upper bound on F k+ i(t), then F k+ i{t) must 
cross it at some point. If this occurs at t = in, then 
F k+1 (t ) = / fc+ i(*o) and F fc+ i(*o) > A+i(*o) 35]. But 
must still satisfy the inequality in equation (|2a|l . 

Thus 

< 2y/F k+ i(ta)y/F k (t ) - F k+1 (t ) 

< 2 V /fc+l(^o) V /fc(*o) - fk+l(to) = fk+l(to), 

which contradicts the assumption that F k+ \ (t) crosses 
f k +\{t) at to- Thus if f k (t) is an upper bound, then so is 
fk+i(t). _ _ 

The initial step in the induction (that fx (t) is an upper 
bound) follows from the second of the rate equations Il2bll , 
and the fact that the generalized singlet fraction is upper 
bounded by 1. Fig.[2shows numerically calculated curves 
fk(t) saturating the rate equations. 




time t 

FIG. 2: Numerically calculated generalized singlet fraction 
curves fk(t) saturating the inequalities in the entanglement 
rate equations 12al2bH . The final solid curve is for k — 50, 
i.e. the singlet fraction of the end qubits in a separated by 
an interaction distance d — 100. (The dashed curves show 
the corresponding upper and lower bounds u k (i) and lk(t) for 
k = 50, from subsection II V Bll . 



IV. HOW FAST CAN ENTANGLEMENT BE 
CREATED? 

How fast can entanglement be generated in a system of 
interacting particles? The question is both theoretically 



interesting, and experimentally important. Many quan- 
tum information processing tasks require entanglement, 
and the faster this can be produced, the less the system 
will suffer from decoherence. Quantum computing al- 
gorithms often generate large amounts of entanglement 
during their execution, so determining how fast entan- 
glement can be generated can also provide bounds on 
algorithm complexity |2lj . 

In this section, we briefly review some existing entan- 
glement generation schemes in the context of the entan- 
glement rate equations derived in the previous section. 
We then investigate what universal limits the rate equa- 
tions put on how fast entanglement can be generated, or 
more precisely, how the time required to entangle two 
particles scales with the size of the system. 



A. Entanglement generation schemes 

How fast entanglement can be generated depends, of 
course, on how we are able to manipulate the system. For 
dcfiniteness, consider entanglement generation in a qubit 
chain. It turns out that measurement is a very power- 
ful resource. If we are able to carry out local operations 
on any qubit, including local measurements and classical 
communication of the outcomes, then the end qubits in 
a chain can be maximally entangled in a time indepen- 
dent of the length of the chain. Though not discussed 
in the context of entanglement generation, Briegel and 
Raussendorf showed that a cluster state can be cre- 
ated in a chain in constant time, and local measurements 
on a cluster state allow a Bell-state to be projected out 
on any desired pair of qubits, including the end pair [23[ . 

The constant scaling assumes we neglect the the time 
required for classical communication of the measurement 
outcomes to the ends of the chain. This can be justified 
on theoretical grounds, since classical communication can 
not create entanglement, and it makes sense to consider 
the interactions as the resource. In many physical imple- 
mentations, it is also reasonable on pragmatic grounds: 
classical communication is usually much easier to imple- 
ment than quantum processes. However, if the interac- 
tions are really the only non-local resource, then classical 
communication must also be implemented via the chain, 
and local measurements arc of no benefit, which is equiv- 
alent to the local-control scenario described below. This 
might be the relevant scenario, for instance, for quantum 
computers. 

If we can apply local unitary operations on any qubit in 
the chain, but not measurements, then we can efficiently 
simulate evolution under any Hamiltonian (this is true 
for general s yste ms of interacting particles, not just for 
qubit chains |24L l25| 'l . Again, it is reasonable to discount 
local resources, which in this scenario means neglecting 
the time required to carry out the local unitaries (the 
'fast local unitary' approximation). And again, this can 
also be justified on physical grounds, since local unitaries 
are typically much faster than interactions. 
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Khaneja and Glaser have developed an interesting pro- 
tocol for state transfer in this scenario [2(| , in the context 
of NMR spectroscopy, which can easily be transformed 
into an entanglement generation protocol. First the mid- 
dle qubits are entangled, then the state of each middle 
qubit is encoded into a three-qubit state. The encoded 
states are transferred along the chain towards the ends, 
where they are decoded again. The protocol requires lo- 
cal unitaries to be applied at discrete times. The evolu- 
tion of the generalized singlet fractions is shown in Fig [31 
clearly reflecting the fact that the protocol is based on 
moving states step-by-step along the chain. It achieves 
a surprising three-fold speedup over the trivial swapping 
protocol for entanglement generation in a chain (entan- 
gle the middle qubits; move to the ends by swapping), 
though the scaling of the time with the length of the 
chain is still linear, as in the trivial protocol. In the next 
subsection, we will use the entanglement rate equations 
to derive a lower bound on the scaling in this local-control 
scenario. 




0.2 



0.1 

„l 1 1 1 1 1 1 1 

01 2345678 

time t / jc/4 

FIG. 3: Entanglement dynamics in the entanglement gener- 
ation protocol based on Khaneja and Glaser's state transfer 
scheme [2f|, for a chain of 10 qubits. Successive curves show 
the evolution of generalized singlet fractions F\ through F5, 
numbered as in subsection 1111 CI 

Finally, we may have no local control over the qubits, 
only retaining the ability to switch on interactions in the 
entire chain, and switch them off at some later time. 
Christandl et. al. developed a state-transfer protocol for 
qubit chains in this scenario [27|. and Yung et. al. have 
given a simple extension to entanglement generation [28| . 
The only local control required is fixing the coupling 
strengths between different qubits, which must be inho- 
mogeneous. Fig. 0] shows the entanglement dynamics for 
the odd chain-length protocol of Ref. [28| — very differ- 
ent to that of Fig. |21 If the strongest coupling strength is 
normalized to some fixed value, then the time to create 
a maximally entangled pair again scales linearly with the 
length of the chain [2!j • 

Osborne and Linden have also developed a protocol for 




0.3 - 
0.2 - 
0.1 - 

Q | 1 1 1 , 

0.5 1 1.5 2 

time t / rc/4 

FIG. 4: Entanglement dynamics in the entanglement gener- 
ation scheme for odd chain lengths from Ref. [28|. here for 
9 qubits. Successive curves show the evolution of general- 
ized singlet fractions F\ through F4, numbered as in subsec- 
tion [^^H] (Note that times are not comparable to those in 
Fig- EH since interaction strengths in |2g| are not normalized). 

state transfer in qubit chains, which could be adapted 
to entanglement generation, involving limited local con- 
trol over a vanishingly small (in the limit of large chain 
lengths) number of qubits at each end of the chain |30| . 

B. Bounds on entanglement generation 

In this subsection, we will use the entanglement rate 
equations to derive a lower bound on how the time to 
create a maximally entangled state scales with the size 
of the system. 

Unfortunately, the set of differential equations defined 
by the rate equations I|2al2b|l has no known closed-form 
analytic solution (at least, none that we could find in the 
literature). Solving numerically can provide numerical 
bounds on the time required for entanglement generation 
(see Fig. [21 . The interesting question, though, is how this 
time scales with the size of the system (for instance, the 
length of a chain), which requires an analytic result. 

For simplicity, we will derive a bound on the scaling of 
the time to entangle the end qubits in a chain of length 
L. The YL/2\ generalized singlet fractions Fk will be 
numbered such that F\l /2j is the singlet fraction of the 
end two qubits. We assume all interaction strengths are 
equal to 1, and that the chain is initially in a completely 
separable pure state (thus Fk(t = 0) = 1/2 for all k). The 
result can easily be generalized to different interaction 
strengths, and indeed to general networks of particles 
(c.f. discussion in subsection IIIIBII . 

We are interested in the time at which F\l/2J (the 
singlet fraction of the two end qubits) reaches 1, as a 
function of L. Though the rate equations do not have 
an analytic solution, we can inductively prove a bound 



8 



on the scaling of this time with L, using an argument 
similar to that used in subsection 1111 CI 

There, we showed that the curves obtained when the 
inequalities in the rate equations are saturated give up- 
per bounds on the evolution of the generalized singlet 
fractions. We can use the same argument to prove that 
we still get upper bounds if we weaken the inequalities, 
using the fact that Fk(t) < 1, and instead solve 



fk(t) = 2VA-1W- AW- 

We can use the argument a third time to prove that if 
Uk(t) is an upper bound on the new fk(t), then the solu- 
tion Uk+l(t) to 



Uk 



+i(t) = 2y / u k (t) - u k+ i(t) 



(3) 



is an upper bound on f k+ i(t). I.e. we have u k {t) > 
fkif) > F k (t). (As concerns boundary conditions, we 
simply require that u fc+ i(0) > /fc+i(0) = F k +i{0) = 
1/2.) 

Now assume there is a u k (t) of the form 



Uk(t) = — 



1 



(4) 



that is an upper bound on f k {t) for some positive con- 
stants a k and e. The differential equation for u k+1 (t) 
then has a solution of the same form as u k (t) (as can 
be seen by direct substitution), with a k +i given by the 
recursion relation 



a k . a k 4 
Ofc+l = — + — a/1 + — . 
2 2 V «fe 

Since Uk+i(0) = (l+e)/2, which is greater than the initial 
condition /fc+i(0) = F k+ i(0) = 1/2, u k+ i(t) is an upper 
bound on f k +i(t) by the argument above. 

All that remains is the initial step in the induction: 
that there is indeed a bound u\{t) on f\(t) with the form 
assumed in Q, for some constants a\ and e. Fortuitously, 
the differential equation (|2b|) for F\(t) (the generalized 
singlet fraction of the entire chain, split into two halves) 
can be solved analytically when the inequality is satu- 
rated (and without weakening the inequality). The solu- 
tion has the form 

/i(t) = sin 2 (* + </>) 

with <fr an arbitrary constant. There is also a trivial so- 
lution: /i(t) = 1. Since the chain starts in a completely 
separable pure state, the initial condition is /i(0) = 1/2, 
and the solution we require is 



fi(t) = 



'sin 2 (t + tt/4) t < tt/4 
1 t > tt/4 



The two parts to the solution merely reflect the fact that 
once the generalized singlet fraction has reached its max- 
imum value of 1 , there is nothing to be gained by further 



interaction, and the interactions affecting F% (namely 
the interactions in the middle of the chain) should be 
switched off. 

Knowing an explicit solution for /i(t), it is easy to find 
a bound u\(t) with the appropriate form. To make the 
algebra simpler, we can upper-bound fi(t) by t + 1/2. 
Thus a ui(t) with the form given in (0J that satisfies 
Ui(t) > t + 1/2 will suffice to complete the proof. This 
leads to the relation a\ < 2e. Any positive a\ and e 
satisfying this will give an appropriate u\{t) > fi(t) > 
Fi(t), and will guarantee that ui(0) = (l+e)/2 > /i(0) = 
-Fi(O) = 1/2. Therefore we have shown that an upper 
bound on F\ (t) with the appropriate form exists, which 
completes the proof. For neatness, we can let e — > 0, so 
that Ufc(0) — > F k (0) = 1/2 and a — > (as used to give 
the curve Uso(i) shown in Fig. [21). 

Solving it[L/2j (t) — 1 gives a lower bound on the time 
required for /[L/2J to reach 1, which is itself a lower 
bound on the time T ent required for the singlet fraction 
of the end two qubits -F\l/2J ^° reacn 1) or equivalently, 
for the end qubits to become maximally entangled. 

We are interested in the scaling of T ent for large chain 
lengths, when a k becomes large. Rather than solving 
U\L/2\(t) — 1 explicitly to obtain the bound, we can 
Taylor expand the square-root in the recursion relation to 
show that it asymptomatically approaches a k = a k -i + 1, 
or equivalently a k = ai + k, as k — > oo. Thus for large 
L, the bound tends to 



To, 



> 



[L/2\ 



a square-root scaling with chain length (see Fig.[5J). 




1 50 200 250 

Chain length L 

FIG. 5: Scaling with chain-length L of the time T ent re- 
quired to create a maximally entangled state between the 
ends. The points show numerical results obtained by saturat- 
ing the rate equations (12al2l>l . The solid and dashed curves 
show the analytic lower and upper bounds, T cnt > \J [L/2\/2 
and T cnt < \J \_L/2\ respectively. 

We have loosened many an inequality during the proof 
of the square-root bound. Could the rate equations give 
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a tighter bound? We can use essentially the same proof 
with the inequalities reversed to prove that a square-root 
bound is the best that can be obtained. 

Instead of using Ff.(t) < 1 to weaken the inequality 
right at the beginning, we use F^{t) > 1/2, which is 
valid when the Fk (t) saturate the inequalities in the rate 
equations, i.e. when Fk(t) = fk{t)- Then, solutions of 

h{t) = y/2y/l k -i(t) - l k (t) 

are lower bounds on fk(t). We can rescale the time r = 
t/y/2 to so that the differential equation for Zfc(r) has the 
same form as that for Uk{t) in the previous proof (Eq.[3J. 
Assuming solutions of the form /^(t) — r 2 /afc, solving 
the resulting recursion relation, and proving there is a 
lower bound on f\(t) of the appropriate form, leads to an 
upper bound on the scaling, for any evolution saturating 
the rate equations. For large chain lengths, the bound 
tends to T ont < yJ\L/2\ — also a square-root scaling. 
Therefore, the square-root bound we have derived is, up 
to a \[2 numerical factor, the best that can be obtained 
from the entanglement rate equations (see Fig. |SJ. 

How does our bound compare with the entanglement 
generation protocols described in the previous subsec- 
tion? The generalized singlet fractions evolve quite dif- 
ferently in those protocols, compared to the evolution 
that would saturate the rate equations (compare Figs. OH 
and 0] with Fig. (2J . All existing protocols that we know 
of scale linearly with the length of the chain — no better 
than the trivial swapping protocol (entangle the middle 
qubits; move to the ends by swapping). It is an inter- 
esting open problem to determine whether any protocol 
can achieve a square-root scaling, or whether the bound 
derived via the rate equations is too weak and can not 
be saturated (which would suggest some improvement on 
the rate equations might be possible). 



V. CONCLUSIONS 

We have investigated entanglement flow, both through 
individual particles and along networks of interacting 
particles. In both cases, we have derived differential 
equations relating the rate of entanglement generation 
to the existing entanglement in the system. 

Entanglement flow through a particle is limited by the 
entanglement of that particle with the rest of the sys- 
tem, providing the system is in a pure state. (Previ- 
ous work |12| has already shown that the entanglement 
can be transmitted by a particle without that particles 
becoming entangled at all, if the system is in a mixed 
state.) 

To describe entanglement flow along general networks 
of interacting particles, we have derived a set of entangle- 
ment rate equations, analogous to the rate equations for 
a chemical reaction. These can intuitively be interpreted 
as describing a flow of entanglement along the network. 
We have used the rate equations to prove a square-root 



lower bound on the scaling with system size of the time 
required to create a maximally entangled state, and com- 
pared this to existing entanglement generation protocols. 
Whether this bound is achievable, or whether the rate 
equations can be improved to give a tighter bound, re- 
mains an interesting open problem. 

The entanglement rate equations were derived in the 
context of two-qubit entanglement creation. However, 
since they involve fidelity-based quantities, they can eas- 
ily be extended to more general settings. Firstly, the 
quantities and equations can be extended to bipartite 
entanglement generation in arbitrary spaces, by taking 
fidelities with a bipartite maximally entangled state in 
the appropriate space. Secondly, they can be generalized 
to the multipartite setting, by taking fidelities with the 
desired multipartite entangled state (e.g. a GHZ state), 
rather than with a bipartite entangled state. 

Together, the results establish a quantitative concept 
of entanglement flow in interacting systems. This is of 
interest as an abstract concept in itself, but could also 
be interesting both theoretically and practically: in the 
analysis of quantum algorithms, for example, since these 
often involve creating large amounts of entanglement dur- 
ing their operation. Or in physical implementations of 
quantum systems, in which it is important to carry out 
any manipulation (including entanglement creation) as 
fast as possible, to beat decoherence. 



APPENDIX A: THREE-QUBIT CHAIN 

To derive the three-qubit result, we use a matrix anal- 
ysis approach to calculating the concurrence, developed 
in [l5j . Writing the Schmidt decomposition of the three- 
qubit system with respect to the partition (b : ac) as 

I^Lc = A i \<Pi)ac + A 2 \<P2) ac |X2> 6 , we can represent 
the state of ab by a 4 x 2 matrix X = (Ai \(fi) , A2 |v?2))- 
The reduced density matrix is then given by p ac = XX' . 

The concurrence C ac of p ac can be obtained from the 
singular values of A = X T T,X, where S = a y ®a y In 
our case, A is a 2 x 2 matrix (because p ac has rank two), 
with just two singular values: > <?2- Thus C ac — ft— ft- 
Since tr A^A — <;f + <rf and | det^4| = ft^, we can also 
write this as 

C 2 ac = trAU-2|detA|. (AI) 

To calculate the time-derivative of the concurrence, we 
must calculate the derivatives of tr(A'A) anc i | det A|. 
From its definition, A = X T Y,X + A T E X. Meanwhile, 
d(tr A^A)/dt = tr(A^A + A^A), which, after a little al- 
gebra, leads to 

d(tr^U) =4Re ^ tr ( S)Q * s ^ x t)). (A2 ) 

Since A is a 2 x 2 matrix, det A = tr(Aa y A T a y )/2. 
Thus d(det A) /At = ti(Aa y A T a y + Aa y A T a y )/2 which, 
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after a little more algebra, gives 



d(det A) 
dt 



Atr[X(T y X T T 1 Xa y X T Y I 



(A3) 



The three-qubit chain evolves according to the Hamil- 
tonian H = H a b g l c + 1„ g Hb c - The two-qubit 
Hamiltonian H a b has a product decomposition H a b — 
Yljj a-ijOi (8 <Jj, where the Pauli matrices Oi are defined 
in the {|xi)Jx2)} basis, and coefficients are real. 
Similarly for H bc and coefficients cy. The Schrodinger 
equation describing the evolution of the system state \tp) 
translates into an equation for the evolution of X: 



X 



ij 



IXai 



c i3 1 ® o-j XaJ) 



We can use this, along with expressions (|A2(I and i|A3|l . 
in the time-derivative of l|Alfl to obtain an expression for 
the derivative of the concurrence: 



r\C 2 



The factor h(H, \i/>)) = J2ij Oijfc*- 



hlj depends on 



ij'Hj 



both the interactions and the system state, and is a rather 
complicated sum over terms involving a%j and cy. We 
define 

t% = ^i\\®a k \vi) 

Oij = {<Pi\<Pj) 

Kx = -i(A?si 2 o u + \ 2 2 s 21 o 22 ) 

h iy = ^2 S 21°22 - ^lS 12 Ou 

h 1z = -iAiAafs^oia - s\ 2 o 21 ), 
and define /i?- similarly to hfj , but with the sy 's replaced 



by ty 's. Note that s 
spin-flip operation [1 



10) 



0. (The tildes denote the 



h(H, ll{e()J -/,,//';, •'•,,/'',,) • 1 ^".J'l •'•,.,//',, 



cr y |</3*).) Then 



However, as we are primarily interested in the de- 
pendence on entanglement (i.e. the dependence on the 
Schmidt coefficients), we can bound the magnitudes of 
the s^-, and Oy by 1, and assume all terms sum in 
phase, giving the bound: 

h(H, |V)) < 8 



E 



\dij\ + |cy| 



which is independent of the system state, depending only 
on the interaction strengths. 



APPENDIX B: ENTANGLEMENT RATE 
EQUATIONS 

The proof of the entanglement rate equations revolves 
around Uhlmann's theorem [HI l3lj . which relates the 



fidelity of two mixed states to the fidelity of their purifi- 
cations: 

Theorem 1 (Uhlmann) If p and a are two states in 
the same Hilbert space Ti, let |V) and \(p) be purifications 
of p and a into a (in general larger) Hilbert space Ti ® Ti' . 
Then 

F(p,a) = max \(<p\ip)f 

where the maximization is over all purifications. 

Since any purification can be transformed into another 
by a unitary acting on Ti', we can fix one of the purifica- 
tions and only maximize over the other one. Also, global 
phases can be chosen to ensure the overlap {(p\ip) is real 
and positive, so the absolute value can be dropped. 

Recall that the entanglement rate equations (|2al2b|) 
involve two disjoint subsets, A and B, of the entire set 
of particles S, which are interacting via two-particle in- 
teractions Hij. We can apply Uhlmann's theorem to the 
generalized singlet fraction of pab at time t: 

FAB{t) = max (0| ti/ ab (U A ® U B Pab (t) U\ ® ul) \<t>) 



Ua,U b 



= max (x| U A ® U B \W by Uhlmann (Bl) 

Ua,Ub 

\x) 

= (x\ Ua ® Ua |^) 2 , 

where \x), Ua and Ub denote the particular state and 
unitaries achieving the maximum. (We are retaining the 
unitaries, rather than incorporating them into one of 
the purifications, for later convenience. Strictly speak- 
ing, they should be extended to Ti ® Ti' and written 
Ua ® Ub <8> Ircst- In the interests of economy, we will 
drop all l res t's.) 

The state Ua <8> Ub |V) ca n be chosen to be any fixed 
purification of the two-qubit density operator a a b- We 
use that freedom to make a purification of the overall 
system state ps, which guarantees that Ua ® Ub \ip) is 
a purification of a a b, as required by Uhlmann's theorem. 
As for \x), since \<j>) is already pure, it is simply an ex- 
tension toTi® Ti': \\) — \4>) ab |$) res t (* ne maximization 
then being over \&)). 

If the system evolves under the Hamiltonian H = 
y~]jj Hij for an infinitesimal time 8t, the state evolves 
to PAB (t + St) = tr /AB (e- iHSt p s (t)e imt ). By writing 
the density matrix of the entire system, ps, in a prod- 
uct basis for the partition (AB : rest) and expanding the 
exponentials to first order in St, it is straightforward to 
show that only interactions involving at least one parti- 
cle in AB give a first-order contribution to the evolution. 
Therefore, H need only include that smaller set of inter- 
actions. Letting 



Ust = exp(-iSty^Hi 



jeAB 
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be the resulting (infinitesimal) unitary evolution opera- 
tor, the singlet fraction after the evolution becomes 

F AB {t + St) = max (C| V A ®V B - U Bt 10" , (B2) 

Va,Vb 

10 

where we have used Uhlmann's theorem again. The state 
\Q is again simply an extension of \<j>) to Ti <8> W, and 
V A <8> Vb Ugt \0 can be chosen to be any fixed purification 
of the two-qubit density operator 

Tab = tr M (V A (8 V B U St ps U ] H Fj ® VjQ , 

Again making use of this freedom, and recalling that we 
chose |V) to be a purification of ps, we can choose |£) to 
be the same state as before: |£) = |V)- 

The state |x) and unitaries L/a and Ub were defined to 
be those maximizing expression IjBlj) . Thus by definition, 

(xl Ox Ub |V) > (xl U A <8> I7 fl |V) 

for all ?7a and C/b. I n particular, this is true for 
infinitesimal changes, e.g. |x) + St Ix' 1 ) where \x^) is or- 
thogonal to Thus (x ± \ Ua <8> C/s |V) < 0. However, 
if this were strictly negative for some Ix' 1 ), then — \x ) 
would make it positive. Therefore (x ± \ Ua ® Ub |V) = 0. 
Similarly, considering infinitesimal changes to the uni- 
taries, we can show that: 

(X| U A H A <8 U B |V) = (xl #A <8 ffuFu |V) = 

(B3) 

(x x l ^ ® % IV) = o. 

Expression l|B2(l for the generalized singlet fraction at 
time t + St must tend to expression (|Bljl (the correspond- 
ing expression for time t) as St — > 0, so |C) = |x) + ^ Ix -1 ) 
and Va (B ) = U A(B ){^ + iStH A(B )), where Ha(b) is a 
Hermitian operator on A{B). Using this, expanding 
Ug t = 1 — iStH + 0(St 2 ) (where H is the sum of in- 
teractions involving at least one particle in A or B), and 
making use of relations (|B3ll , we have 

F AB (t + St) = (x| U A <g> U B (1 - iStH) |V) 2 + O (St 2 ) . 

I.e. the state and unitaries maximizing expression i|Bl|) 
also maximize (IB2() . to first order in St. 

Hamiltonian H currently includes all interactions in- 
volving at least one particle in ^4 or B. By expanding 
H in the previous expression as a sum over these two- 
particle interactions, we can use the same relations <|B3|) 
to show that only interactions crossing the boundary of 
A or B need to be included to give the generalized singlet 
fraction to first order in St: 

F AB (t+St) = (x| U A ®U B (l-i6tJ2 H ij) Hf+O (St 2 ) . 

i£A(B) 
j£A(B) 

Now, global phases were chos en to m ake it real and 
positive, so (x| U A ® Ub |V) — \fF AB {t). Thus, expand- 
ing the square in the previous expression and only retain- 
ing first order terms in St, we arrive at a first expression 



for the time-derivative of the generalized singlet fraction: 
F AB (t) = ^FMtj- X - (<p\ H tj |V) - (VI H tj \<p) (B4) 

i£A{B) 
jeA(B) 

where \<p) = u\®U B \ X ). 

To proceed, we will need the following Proposition |33| 
which we use to prove the subsequent Lemma: 

Proposition 2 (Fan-Hoffman) For any operator X , 
the ordered singular values o~\ of X are individually 
greater than or equal to the ordered eigenvalues r\ of 
ReX = (X + Xt)/2. I.e. a\ > r\ V*. 

Note that the eigenvalues of ReX can be negative, in 
which case the absolute values of the eigenvalues need 
not obey the Proposition. 

Lemma 3 For any operator X , (tr |A|) 2 — (tr(ReA)) 2 > 
tr((ImA) 2 ), where \X\ = y/xlO, ReX = (X + X r )/2, 
and lmX = (X - X*)/2i. 

Proof Assume initially that tr(ReX) is non-negative. 
Defining P (N) to be the set of positive (negative) eigen- 
values of Re X , 

= ^a i a j - r * r i + E W W 

+ ENN - E MN 

=E ff ^- + ENN+ E '•;':;) 

iorj^P i£P ijeP 

+en(en-ei4 

ieN S'eP j&N ' 

The first two terms are clearly positive, the third is pos- 
itive by Proposition [21 and the last by the assumption 
that tr(ReA) > 0. Thus £\ (a^j - nrj) > 0. Now, 

(tr|A|) 2 - (tr(ReA)) 2 = a t ) " - (E n) * 

i i 

= E a * ~ E r ^ + 2 E (? ia i ~ 
>E^-E-? 

i i 

= tr(AA t ) -tr((ReA) 2 ) =tr((ImX) 2 ) 

where in the last line we have expanded X = Re A + 
ilmA, and used the fact that Re A and ImA are both 
Hermitian and that the trace of their commutator is 0. 
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For completeness, we can remove the assumption 
tr(ReX) > by noting that, if there existed an operator 
X with tr(ReX) < such that the Lemma did not hold, 
then the operator —X would also violate the Lemma. 
But then tr(Re(— X)) > 0, so the Lemma must hold for 
all operators. □ 

Recall that Hij really means iJy ® l r cst- Thus 

Ha |V) - h.c.) = tv^Hij ■ j (tr Aj |V)<^| - h.c.)) 
= 2tr {Hij YvaXij) where Xij = tvuj \ip)(ip\ 



< 2,/tr Hf, Jtr((Im ) 2 ) by Cauchy-Schwartz 



< 2 II H„ 



y ii hs 



.y(tr|^|)2-(tr(Re^-)) 5 



(B5) 



using Lemma |3 in the last line. (|| • ||hs denotes the 
Hilbert-Schmidt norm.) 

Finally, we need to relate the quantities under the 
square-root to generalized singlet fractions. Firstly, 
(tr(ReA%)) 2 = (Rc(tr|^|)) 2 = {^f = F AB (t), 
since global phases were chosen to make ((p\tp) real and 
positive. 

Secondly, acts on one particle j within A or B, and 
a particle i outside. If j is in A, define the sets A[ = AUi 
and B[ = B. If it is in B, define A\ = A, B[ = B U i. 
We apply Uhlmann's theorem to the definition of the 
generalized singlet fraction for pa'. B '., an d again choose 
the same state for one of the purifications. So long 
as A\ and B[ are disjoint, we have 



F 



A[B[ 



(i) = max 

v A ,,v B 
10 



(CI V A > ® Vb' 



> max 

v.,,v„, 



(X\ U A V A > ® UbV b > IV') 2 
(<p\ V* (8 Vb' |^) 2 > max((^| tf y |^) : 



= max 
v A ,,v E 

= m r ax(tr(^ i Xy)) 2 =tr(|Xy|) 



where an inequality appears each time we restrict the 
maximization. The last line follows from the fact that, 
for any operator, tr \X\ = maxj/ | tr(UX)\, which is easily 
proved via the polar decomposition of X. If A\ and B\ 
have a particle in common (it must be particle i if they 
do), then the second of the two inequalities is not valid. 
We can instead bound tr(|A i:) -|) 2 < 1. 

Thus, using (tr(ReX y )) 2 = F AB (t) and tr(|Xy|) 2 < 
F A[B[{t) {A and B[ disjoint) or tr(|X^|) 2 < 1 {A\ and B[ 
overlapping) in (|B5|) . and substituting the result in (|B4|) . 
we arrive at a version of the entanglement rate equation: 



Fab (t) < 2 



i£A(B) 

jeA(B) 



F AB {t)JF A , B ,{t)-F AB {t), 



where F A i B i is defined to be 1 if A\ and B[ overlap. 

To describe entanglement flow in a network of inter- 
acting particles, one such rate equation must be written 



down for all meaningful generalized singlet fractions that 
can be defined on the network ('meaningful' implying 
that the sets A and B include particles a and b respec- 
tively, and are each be made up of 'one piece'). 

Recall that, since A\ and B[ are subsets of A' and B' 
(see Fig. QJ, F A > B > < F A / B /, We can use this to arrive 
at the simpler version of the entanglement rate equations 
presented in the main text: 



Fab® < VF AB (tWF A , B ,(t) - F AB (t). 

i<£A(B) 
jeA(B) 



APPENDIX C: GENERAL TRIPARTITE CHAINS 

The first half of the derivation of the entanglement 
rate equations given in Appendix [5] can be re-used in 
the proof of the tripartite chain result. Recall that the 
three systems A, B and C making up the chain are in 
an overall pure state |V)abC' and interact by nearest- 
neighbour interactions: H = H AB + H B c- As noted in 
subsection IIII Al the entangled fraction F A c of p A c can 
be expressed as a maximization over unitaries U A and 
Uc rather than states. Applying Uhlmann's relation, it 
can be rewritten 

F AC = max (x| U A ® U c \ip) = (</#) ■ 

Ua,Uc 

\x) 

We can choose \ip) to be the overall system state |V)j4bc> 
since this is a purification of pac- \x) i s then an ex- 
tension of a maximally entangled state \<j)) A( j on 'H.ac 
to the space Habc- \x) = \4>)ac I Ob- We define W) = 
)ac I0> where bars denote the particular 



U 



t U-A - 



w c \x) 



unitaries and states achieving the maximum. 

Although the rate equations involve the generalized 
singlet fraction, up to expression (|B4I) the derivation in 
Appendix [H] applies equally well to the entangled frac- 
tion. Expression (|B4(1 then becomes 



F AC (t) 



y/F^(fi.l({<p\H\1>)-{1>\H\<p)), (CI) 



Now, writing the state in its Schmidt decomposition 



'/AC l"/B> 



for the partition (B : AC), \'4>)abC 
where we sort the Schmidt coefficients in descending or- 
der: Ai > A2 > • ■ ■ > A n . Extending {|«) B } to form a 
complete basis for Tl B , \ip) can be written in the product 

decomposition \ip) = \4>)ac Ei a i N)b (^ ne ai are com pl ex 
in general). 

We know that \(p) maximizes (ip\ip) = Ej a i ^» (^IVi)- 
Clearly, the phases of a, must be chosen to cancel the 
phases of (cf)\ipi). The relationship between the mag- 
nitudes of the ai's and (0|Vi)'s can be found using 
Lagrange multipliers, with the normalization constraint 
Ei l"i| 2 = 1> yielding 



Ai (ipi\4>) 



E fc A 2 . \(4>\t^k)\ 



(C2) 
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We also know that, for any Hamiltonian acting only on 
A or C, (<p\ Ha \tp) = (<p\ H c \ip) — (applying the same 
reasoning as used to prove relations l|B3|l in Appendix[5J| • 
Thus 

«JAi<0| H Ma) \i>i) = -J2 < X ^\ H MC) \1>i) • (C3) 
ijti 

Now, the system Hamiltonian H = Hab + Hac- 
(<p\ H AB IV) = E« <h $1 \j)b l^>- For the 

i = j = 1 terms in the sum, b (1|JJab|1) b is just 
some Hamiltonian acting only on A. Similarly for 
B (l\H B c Thus usin S and 

fo>| # |V) = ^ ajAj ((0, t| ff |Vi, i) - 1| H |^, 1)) 
+ ^2a*X i (4>,i\H\ip j ,j) 
= S ' ft. j j + S ' Aj Aj /ijj , 



where 

h _ (^|0)((0, t |g|^,i)-(^,l|g|^,l)) 

_ {i>S){^i\H\^j) 
<Hj — I = = tor 1 T ]■ 

Using this in (|C1|) . and bounding (hij—h*j)/i < 2\H\ = 
2(\H A b\oo + \Hbc\oo) (where \M\oo = maxy |My| de- 
notes the norm), we arrive at the final result: 

Fac(<) < 2 |ff | V^UcW (Z! X * X i - A i)- 
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